Lande-like formula for the g factors of hole-nanowire subband edges 
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We have analyzed theoretically the Zeeman splitting of hole-quantum-wire subband edges. As is typical for 
any bound state, their g factor depends on both an intrinsic g factor of the material and an additional contribution 
arising from a finite bound-state orbital angular momentum. We discuss the quantum-confinement-induced 
interplay between bulk-material and orbital effects, which is nontrivial due to the presence of strong spin-orbit 
coupling. A compact analytical formula is provided that elucidates this interplay and can be useful for predicting 
Zeeman splitting in generic hole-wire geometries. 
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The interaction of a particle's magnetic moment with an 
external magnetic field B typically results in the lifting of de- 
generacies in its energy spectrum. The study of such Zeeman 
splittings provides important clues about microscopic proper- 
ties of a quantum system. ' In vacuum or the bulk of a solid, 
the Zeeman splitting can be understood as arising from the 
coupling of a (quasi-)free particle's intrinsic angular momen- 
tum (i.e., spin) to B. However, in general, orbital motion 
in a quantum bound state is associated with a magnetic mo- 
ment that interacts with the magnetic field and, thus, also con- 
tributes to Zeeman splitting and the g factor. In the presence of 
spin-orbit coupling, the interplay between orbital and intrinsic 
(spin) contributions to g can be nontrivial.— 

Recently, p-type semiconductor nanostructures have be- 
come available as interesting laboratories to explore the in- 
terplay between quantum confinement and the magnetic mo- 
ment of spin-3/2 holes. ^ ••^•^•^'^'^■^■^ For example, strong spin- 
orbit coupling present in the upper-most valence band of typ- 
ical semiconductors results in an energy splitting'*' between 
heavy-hole (HH) and light-hole (LH) quantum-well subband 
edges." This HH-LH splitting turns out to freeze the HH 
spin-projection axis to be perpendicular to the quantum well, 
thus preventing its reorientation by an in-plane magnetic field 
and, therefore, suppressing Zeeman splitting. ^ In addition 
to HH-LH splitting, confinement typically induces a mixing 
of HH and LH states in quantum wires, '■'■''*''^''^ even at sub- 
band edges. Both HH-LH mixing and splitting strongly affect 
the physical properties of hole states. 

We have previously shown'^ '^ how HH-LH mixing mod- 
ifies the Zeeman splitting of hole-wire subband-edges when 
the bulk-material hole g factor n is large and can be assumed 
to dominate over all orbital magnetic-field effects."^ Here we 
complement this analysis of Zeeman splitting in hole quan- 
tum wires, focusing on the confinement-induced interplay be- 
tween bulk-material and orbital bound-state contributions to 
the g factor. Bulk-material and orbital contributions turn out to 
vary strongly between different quantum-wire subband edges, 
in both their magnitude and sign. Depending on their rela- 
tive sign, the two contributions can either enhance or suppress 
each other, thus yielding large variations in sign and magni- 
tude of the total g factor We provide a formula from which 
the origins of the observed spin-splitting characteristics can 
be explained and which may also be useful to predict Zee- 



man splittings in generic hole-wire geometries. For cylin- 
drical hole wires, we also find a Lande-like formula, which 
explicitly expresses the subband-edge g factor in terms of a 
new total angular momentum and spin-3/2 tensor invariants,— 
the latter resulting from a multipole expansion of the spin- 
3/2 density matrix. It is an interesting feature of spin-3/2 
systems that, in addition to a monopole and dipole moment, 
which are related to the density and polarization, respectively, 
a quadrupole and an octupole moment can exist, which have 
no equivalents in spin- 1/2 systems. Orbital contributions turn 
out to give rise to a dependence of the hole-wire g factors on 
spin-3/2 quadrupole and octupole moments. 

Figure [T] shows our calculated absolute values of buUc- 



wire subband edges. The integer a enumerates wire subbands. 



material and orbital contributions [.g^ I' \9mb \ '■^^ X.oX.?i\ g- 
factor \g[oi\ = |(72"''+5orb^ I' respectively, for cylindrical-hole- 
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FIG. 1: (Color online) Absolute values of orbital and bulk-material 
Zeeman contributions, g'^^"^ and g'^K to the total g-factor, g[^l^ = 
g^"' + ffoit' , for a cylindrical hole nanowire defined by a hard-wall 
confinement in GaAs (i.e., a material with Luttinger parameters 71 = 
6.98, 7s = 2.58, and k = 1.2; see text). Values are shown for 
the highest-in-energy subband edges enumerated by q. fa is the 
quantum number of the total angular momentum operator Fz = ,Jz + 
Lz- The magnetic field is assumed to be parallel to the wire {z) axis. 
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TABLE I: g-factor contributions for cylindrical-hole-nanowire subband edges in GaAs, witii cubic corrections neglected. 
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starting with the highest. The magnetic field is assumed to 
be appHed parallel to the wire, which is aligned with the z- 
coordinate axis, and we used materials parameters character- 
istic for GaAs^*^ but with band warping neglected. Several 
observations can be made, i) The orbital contributions gener- 
ally have a magnitude of the same order as that of the corre- 
sponding bulk-material ones. In some instances, such as for 
levels a = 5, 6 and 10, they greatly exceed the bulk- material 
contributions, giving rise to very large spin splittings, ii) The 
bulk-material and orbital contributions fluctuate in sign and 
magnitude between different subbands. iii) The two contribu- 
tions can variably enhance or suppress each other, thus render- 
ing large variations in sign and magnitude of the total g-factor 

gto" as a function of subband index a. Further below we dis- 
cuss how the observed behavior results from HH-LH mixing. 

Before describing details of our calculational method and 
giving further results, we discuss the general form of the var- 
ious g-factor contributions. As the hole-wire subband edges 
are well-separated in energy, it is possible to calculate their 
Zeeman splitting analytically within a perturbative approach 
that becomes exact in the limit of B = B^z 0. We find 
the following general expression = .g^"^ + 5^,^^diag + 

5orb^mix + forb^cub' '^^^ ^^^^ three tcrms are independent of the 
wire's direction with respect to crystallographic axes, neglect 
band warping, and have the following form: 
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Here (. . . )« denotes the expectation value of the operator in- 
side the angular brackets taken for the positive-helicity eigen- 
state'"*'^' of wire subband edge a obtained in zero magnetic 
field. K is the bulk-hole g-factor^S = xky — yk^, x± = 
X ± iy, k± — kx i: iky, J± = {Jx ± iJy)/^/2, ki are or- 
bital linear momenta, Ji are spin-3/2 angular-momentum op- 
erators, and 7s — (272 + 373)/5 in terms of Luttinger param- 
etersi22i21 The band-warping correction gl"^j^.^^^, depends on the 
wire direction. For now, we continue neglecting this term in 
our calculations, postponing a discussion of its effect on our 
results to the end of the paper. 

As previously discussed, it is a general characteristic of 
hole nanowire subband edges that they are mixtures of HH 
and LH characters i'^'''*i'^ Thus, all of the individual terms in 
the expression of ^to't'' affected by HH-LH mixing, be- 



cause the expectation value is taken with respect to a HH-LH 
mixed nanowire subband-edge state. It was previously shown 
that the bulk-material contribution g^"'' fluctuates strongly be- 
tween different subband edges as a result thereofji^ 

The orbital contributions g^'^fj^i.^^ and ffo^^^ix show similar 
behaviour, as can be seen from their values given in Table U 
Interestingly, with only few exceptions (e.g., the level with 
a — 6), these two terms appear to counteract each other, 
sometimes quite dramatically. (See e.g., levels with a — 4 and 
8.) .goi^'jjiag has a straightforward interpretation in that it rep- 
resents the naive coupling of a bound-state orbital magnetic 
moment to the external magnetic field, with different gyro- 
magnetic ratio for HH and LH states because of their different 

effective masses. The term gl°^^i^ arises because the spin- 
orbit coupling in the valence band induces linear-in-i? terms 
that couple HH and LH amplitudes directly. Thus it represents 
a HH-LH mixing effect on top of the mixing already present 
in the spin-3/2 subband-edge states. 

The total g factor g^^"^ is a sum of bulk-material and orbital 
terms. As seen in Table |I] these different terms can vary sub- 
stantially in magnitude and sign between different subband 
edges a as a result of the HH-LH mixing. More importantly, 
the relative sign between the various contributions is differ- 
ent between the subband edges. Thus, the bulk-material and 
orbital terms can varyingly enhance or suppress each other's 
contribution to the total g factor depending on the subband 
index a. 

We now proceed to discussing in detail the theoretical for- 
malism for calculating the total g-factor and its individual con- 
tributions given in Eqs. (fTa]i-(fTcli. Subsequently, we will show 
that their basic characteristics discussed above are universal 
for quasi- ID hole systems by considering the influence of re- 
duced symmetry of the confining potential. We will also ana- 
lyze the effect of (so far neglected) cubic corrections. 

Our starting point is the Luttinger Hamiltonian^^ in the 
spherical approximation^^ 



Hi^^-—k 14x4+ 7s 



k.J)'--pl 

4 



4x4 



(2) 



In Eq. dU, we used atomic units. The spin-orbit-coupling term 
proportional to 7s is the origin of HH-LH splitting and mixing, 
giving rise to the peculiar spin properties of confined holes. In 
the following, we assume k^ — because we will focus solely 
on the properties of quasi-lD subband edges. The bulk-hole 
Zeeman interaction with a magnetic field applied along the 
wire (z) direction is given by Hz = —2kiibBzJz- (We ne- 
glect the small anisotropic part'^ of the Zeeman splitting in 
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the bulk valence band.) This term is the origin of the contri- 
bution 52 to the hole-wire g factor. The orbital effect of the 
magnetic field is included by replacing k — > k + A, with the 
symmetric gauge A = {—y/2,x/2,0)Bz- As we will only 
consider the leading (linear-in-i?^) contribution to the Zeeman 



splitting, we neglect orbital terms of higher than linear order 
in Bz- The additional Hamiltonian, i/orb^ which incorporates 
the effects due to the interaction between the holes' orbital 
degrees of freedom and the magnetic field, is then given by 



J 



Hoib — —f^sBz 



7l]l4x4 + 7s 



Jz - ^l4x4 



+ ijs { X-k^j'^ — x+k+Jt 



(3) 



r 



The first term, proportional to Lz and diagonal in spin-3/2 



space, gives rise to the contribution gl°^^i^„- 

(a) 

couples HH and LH states and leads to .goit, Tnix- 

For kz = 0, the Hamiltonian Hb = H]^ + Hz + -ffoib 
turns out to be block-diagonal in spin-3/2 space such that 
subspaces for helicity ± do not couple.''*-^' Subband ener- 
gies Ea± within the ± subspaces are found by diagonalising 
corresponding blocks in Hb + V{x, y), where V{x, y) is the 
confining potential defining the wire. For any given a, eigen- 
states |a±) with energies Ea± are found to be degenerate at 
zero magnetic field, but these split for Bz > 0. We define the 
g factor for a hole-wire subband edge with index a accord- 
ing to2^ 5tm' = \imB^^o[Ea+{Bz) - Ea-{Bz)])/{HBBz). 

Because of the finite energy separation of subband edges for 
different a, this definition turns out to be equivalent to the 
perturbative (to first order in Bz) result 



The second term 



(a) 

5tot 



(Hj. + H 



orb / (1+ 



(Hy. + H„, 



2(i/z 



^.■rBz 

Hoib)a+ 



^.■rBz 
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where the expectation values are taken with respect to the wire 
eigenstates that diagonalize the zero-field, quantum-confined 
Luttinger Hamiltonian H]^ + y (x, j/). 

An interesting case arises for hole wires with cylindrical 
cross-section. There, subband-edge bound states are also 
eigenstates of a new total angular momentum Fz = Jz + 
Lz This makes it possible to rewrite Eqs. ([1]) in terms of 

the quantum number /„ associated with Fz, as well as expec- 
tation values of spin-3/2 operators, thereby essentially elimi- 
nating any dependence on the orbital angular momentum Lz- 



In the spherical approximation, we find g^ 

(a) (a) 
?q,cyl + .9o,cyl 



.9r,cyi' where 
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Jz 
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(5a) 
(5b) 

(5c) 
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Here r is the radial coordinate in the xy plane, and expecta- 
tion values are to be taken in the radial spinor wave functions 
that are eigenfunctions of the canonically transformed Hamil- 
tonian e''^-f[Hi^ + V{r)]e-'-'-'^M 

In Eqs. (|5^-c), expectation values are taken of spherical ten- 
sor operators associated with multipole invariants of the spin- 
3/2 density matrix.'^ The contribution .gd.cyi depends on the 
dipole moment, or polarization, of the hole-wire bound state. 
It contains the previously analyzedi^iii bulk-material contri- 
bution that is proportional to k. Interestingly, expectation 
values of tensor components for spin-3/2 quadrupole and oc- 
tupole moments appear in the additional terms ^q.cyi and go.cyb 
respectively. Such quadrupole and octupole moments are a 
unique feature of spin-3/2 systems and have no equivalent in 
the spin- 1/2 realm. Situations can arise in which the dipole 
moment (polarization) is very small, but higher moments can 
be substantial. For example, a magnetic field applied parallel 
to a two-dimensional hole system can induce a large octupole 
moment, whilst the dipole moment associated with the polar- 
ization is vanishing.'^ For the hole nanowires, the appearance 
of higher moments in the g factor are entirely due to orbital ef- 
fects, which is signified by their dependence on the Luttinger 
parameters. The remainder term gi- ^yi is of purely orbital ori- 
gin but does not appear to have a straightforward relation to 
pure spin-3/2 tensors. 

To verify the validity of our approximations and the use 
of the perturbation approach, we have performed calculations 
taking into account the effects of lower symmetries of the con- 
fining potential, as well as the underlying crystal, and com- 
pared these results with numerical diagonalization of the full 
Luttinger Hamiltonian including band warping22i22i2i Our re- 
sults are shown in Fig.|2l where we plot the total g-factors for 
the ten highest wire subband edges with index a = 1 ... 10, 
calculated using a hierarchy of four different approximations: 
I Perturbation theory for iJorb + Hz using unperturbed states 
that diagonalize the Hamiltonian Hl + V{x, y) for a cylindri- 
cal hard-wall potential. (These are the same results as shown 
in Fig. [T]and Table II Same as I, except that the wire is 
defined by a hard-wall potential with square crossection, thus 
reducing the symmetry of the confining potential. Ill Same 
as II, but the unperturbed states are obtained by diagonalizing 
the zreo-field Luttinger Hamiltonian including band-warping 
terms, thus taking into account the lower (cubic) crystal sym- 
metry.^^ (For this case, the wire is assumed to be parallel to 
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FIG. 2: (Color online) Total g-factors for hole wire subband edges la- 
belled by index a = 1 ... 10, calculated using a hierarchy of approx- 
imations described in the text. A close agreement between the differ- 
ent results for low-index subbands is observed, while corresponding 
numbers show an increasing spread for higher-index subband edges. 



the [001] direction.) IV Results obtained from full numer- 
ical diagonalization of the Luttinger Hamiltonian with finite 
magnetic field, including band warping, for a square cross- 
section wire parallel to the [001] direction. ^"^ Comparison be- 
tween these different cases shows good agreement for states 



with low subband index a. For higher a, however, the effects 
of reduced symmetry are more pronounced, and there appears 
an increased spread between the ^-factors calculated for dif- 
ferent approximations. We expect that these will become even 
more important for wires aligned with crystallographic direc- 
tions other than the high-symmetry [001] direction considered 
in our calculations. 

In conclusion, we have presented a comprehensive and 
general^*^ formalism for the description of the Zeeman split- 
ting of hole-quantum-wire subband edges and applied it to 
hole wires subject to a magnetic field applied parallel to the 
wire axis. We have elucidated the role of HH-LH mixing 
in hole nanowire systems, as well as the intricate quantum- 
confinement-induced interplay between the bulk-material and 
orbital contributions to the total g factors. The latter may be 
used for tailoring of hole spin splittings by tuning the rela- 
tive importance of the two contributions, e.g., by electrostatic 
confinement engineering in p-type quantum point contacts^ 
or nanowiresi^ 
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